Spatial noise correlations of a chain of ultracold fermions - A numerical study 
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We present a numerical study of noise correlations, i.e., density-density correlations in momentum 
space, in the extended fermionic Hubbard model in one dimension. In experiments with ultracold 
atoms, these noise correlations can be extracted from time-of- flight images of the expanding cloud. 
Using the density-matrix renormalization group method to investigate the Hubbard model at various 
fillings and interactions, we confirm that the noise correlations contain full information on the most 
important fluctuations present in the system. We point out the importance of the sum rules fulfilled 
by the noise correlations and show that they yield nonsingular structures beyond the predictions of 
bosonization approaches. Noise correlations can thus serve as a universal probe of order and can be 
used to characterize the many-body states of cold atoms in optical lattices. 
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I. INTRODUCTION 

Recent advances in methods for trapping and control- 
ling ultracold atoms have opened up the promising pos- 
sibility of directly simulating strongly interacting many- 
body Hamiltonians. These experimental systems can be 
used to engineer and analyze models that lie beyond the 
scope of present analytical and numerical methods, thus 
potentially shedding new light on fundamental quantum 
many-body problems [1, 2, 3, 4, 5] such as the the na- 
ture of the mechanism for high-Tc superconductivity or 
whether a true spin liquid can be realized. Such systems 
can be formed by trapping atomic gases in one-, two-, or 
three-dimensional optical lattices. The interactions can 
be accurately tuned by adjusting external fields. In a 
similar spirit, it has recently been proposed to simulate 
strongly correlated systems in experiments by studying 
the dynamics of polaritons in arrays of electromagnetic 
cavities, see Ref. 6. 

The techniques for manipulating ultracold atoms are 
fairly advanced and have already enabled a broad range 
of astonishing systems, such as a superfiuid, a Mott in- 
sulator [7, 8] , a strongly interacting Fermi gas [9] or also 
mixtures of bosonic and fermionic gases [10]. However, 
the subsequent analysis of their properties has proven 
to be diflftcult. In view of the application to solid state 
problems, it is crucial to have tools at hand that can 
accurately describe the engineered state, preferably by 
extracting the correlation functions of the atomic gas. A 
recent proposal [11] for a universal probe of correlations 
suggested measuring the shot noise in time-of-fiight im- 
ages of the expanding cloud of atoms after release from 
the trap. This method is based on the fact that, after 
a long enough time of fiight the density distribution 
of the expanding cloud becomes proportional to the mo- 
mentum distribution in the interacting system [11, 12], 



with momentum q (r) = mr/ {ht) for an atom of mass m. 
The noise in the image-by-image statistics is governed by 
higher order correlations of the initial state 

Gaa' (r,r') = (nq(^)nq(^/)) - {n^{r)) (^q(rO> ^ (l) 

where a is an internal quantum number, e.g., the spin, 
that allows different states to be distinguished. By ana- 
lyzing the shot noise in several mean-field states, Altman 
et al. [11] showed that the presence of a particular order 
leaves a very distinctive fingerprint on the noise correla- 
tions, e.g., due to superconductivity or spin order. On the 
experimental side, this quantity has already been mea- 
sured on several occasions, in both fermionic and bosonic 
cold atomic gases, i.e., in bosonic Mott insulators [13, 14], 
fermionic superfiuids [15], and band insulators [16]. 

In the following, we will treat one-dimensional (ID) 
fermionic systems on an optical lattice and will concen- 
trate on the noise correlations of the lattice model itself: 

G{k,k') = Y. Gaa'{k,k') . (2) 

Here the brackets denote the ground state expectation 
value. Shortly after the pioneering analysis of Ref. 11, 
which was based on mean-field calculations, Mathey et 
al. [17] analyzed the shot noise for a ID Tomonaga- 
Luttinger (TL) liquid within a bosonization approach, 
allowing them to explore the momentum-dependence 
around opposite Fermi points, i.e., around k ^ kp and 
k' ^ —kp' In these TL liquids, different types of order 
compete, leading to a rich structure in the noise corre- 
lations. In the present paper, we employ the density- 
matrix renormalization group (DMRG) method [18, 19] 
to study the noise correlations in the ID extended Hub- 
bard model. This numerical approach allows us to go 
beyond the Luttinger theory and to uncover the full set 
of features that are contained in the noise correlations 
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within the entire Brihouin zone. In the vicinity of op- 
posite Fermi points, we find perfect agreement with the 
analytical predictions of Ref. 17. 

The remainder of the paper is organized as follows: 
In Sec. II, we discuss general properties of noise corre- 
lations, independent of the microscopic model. The ex- 
tended Hubbard model considered in this work is then 
introduced in Sec. Ill, together with a summary of the 
different (quasi-) orders encountered in our numerical ap- 
proach. Our main results, the analysis of the noise cor- 
relations for different phases of the extended Hubbard 
model, are presented in Sec. IV. In Sec. V, we make the 
connection to atomic physics and discuss experimental 
issues. We present our conclusions in Sec. VI. 



II. PROPERTIES OF NOISE CORRELATIONS 
IN FERMIONIC LATTICE MODELS 

Before considering a specific microscopic model, it is 
useful to discuss several general properties of the noise 
correlation functions at the lattice level. In the follow- 
ing, we consider the noise correlations (2) on a periodic 
lattice with L sites. In this case, k denotes the lattice 
momentum and a describes an internal quantum num- 
ber that we associate with the spin, i.e., a G {T? i}- 
general, a can denote a more general fiavor or species 
index; the statements below remain valid as long as the 
density operator 



can be written as a product of creation and annihilation 
operators that satisfy the canonical fermionic commuta- 
tion rules 

It is convenient to use the Fourier transformation 



(3) 



SO that the creation and annihilation operators in coordi- 
nate space also obey the standard commutation relations. 
Under these assumptions, G satisfies the following exact 
statements: 

1. Bounds: For all k^k' and cr, a', G is uniformly 
bounded 

1 



\G^a' {k,k')\ < 



(4) 



independent of the system size, i.e., Gaa'{k^k^) it- 
self cannot diverge with system size. 

2. Sum rules: In a system in which the number of par- 
ticles of every species is conserved, the noise corre- 
lations satisfy the sum rule 

{k,k') = {N^N,.) - (N^) (N,.) = , (5) 

k,k' 



where = J2k ^k,a- 

3. Equal-spin momentum diagonal: Along the mo- 
mentum diagonal k = k' ^ the equal spin cj = cj' 
shot noise is given by 



G,a(/c,/c) = (nfc,,)-(nfc,,)2>0. 



(6) 



This part of the noise correlation function is there- 
fore completely determined by the momentum dis- 
tribution function (uk^a) itself and does not contain 
any other information. 

The basic idea behind the use of density-density cor- 
relations as a universal probe of the system lies in the 
simple fact that they contain important parts of various 
particle-hole and particle-particle correlation functions. 
Quite generally, a particle-hole scattering operator with 
net momentum q and form factor f{k) can be written as 



(7) 



with 



k,a,f3 



for generalized charge-density waves (CDW) and spin- 
density waves (SDW), respectively [20]. Here a are the 
Pauli matrices. This ensemble of CDW and SDW op- 
erators encompasses a large variety of particle-hole in- 
stabilities, such as conventional charge and spin-density 
waves, bond-order waves, staggered fiux phases, as well 
as spin nematic states [20]. In ID fermionic systems, the 
strongest instabilities in the particle-hole sector often oc- 
cur at q = 2kF' However, in some cases, the g ^ 
correlations can also become strong, for example, close 
to a ferromagnetic transition, or when longer-range in- 
teractions induce important correlations at wave vectors 
different from 2/ci?. Analogously, in the particle-particle 
sector, a scattering operator of particle pairs with total 
momentum q and form factor / reads 



1 



^ k,a,f3 



f{k)af3 cl 



k-\-q,a'-k,l3 



(8) 



Here the dependence of the form factor on the spin 
indices determines whether the operator describes 
spin-singlet (SS) or spin-triplet (TS) superconducting 
pairing. In the case of SS, / is an even function in the 
momentum /c, whereas for TS, / is odd in k. In the fol- 
lowing, we restrict our discussion to the most common 
case of pairs with total momentum (7 = 0, but we remark 
that the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) [21] 
state appearing in the current discussion on the super- 
fiuid pairing properties in population-imbalanced Fermi 
gases would require a finite q. 

The strongest correlations can be identified by compar- 
ing the structure factors of possible competing orderings. 
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TABLE I: Dominant correlations for the parameter sets studied in this work, supplemented with schematic pictures of the 
ground states. LRO stands for true long-range order, while QLRO denotes algebraic correlations (quasi-long-range order). 
We encounter charge- density waves (CDW), spin-density waves (SDW), bond-order waves (pCDW, where the kinetic energy is 
modulated), singlet superconductivity (SS) and triplet superconductivity (TS). The brackets indicate the subdominant correlation 
functions, i.e., ones which decay with a faster power law than the dominant correlations. The momentum of the particle-hole 
correlations is 2kF — 27tx, unless stated otherwise. The pairing correlations are for pairs with zero total momentum. 



With the above definitions of the scattering operators, 
the associated structure factors can be expressed as 



Si{f,q) = (oUf,q)Oi{f,q) 



(9) 



In the case of true long-range order, the structure fac- 
tor taken at the ordering wave vector k would grow lin- 
early with the system size and would thus diverge in the 
thermodynamic limit. For the quasi- long-range order en- 
countered in many ID systems, in which the asymptotic 
decay of the correlation functions has the form of a power 
law (algebraic decay) , the corresponding structure factors 
can show similar power-law divergences, with exponents 
that depend on the value of the interactions in the micro- 
scopic model. To see the tight correspondence between 
the structure factor (9) and the noise correlation func- 
tion (2), it is illuminating to rewrite the operators in a 
slightly different way. For instance, the four-body opera- 
tors of the particle-hole structure factor that contribute 
to the noise correlations read 



/(fc)a/?/(fc')a'/?'4 



k'=k 
a=a' 

f3=f3' 

2 „t t 



k-\-q,a'^k,/3'^k' ,/3''^k'-\-q,a' 



where the expression on the second line is easily recog- 
nized as one term of the noise correlations. The addi- 
tional terms are two-body operators that arise due to 
commutation rules. From this "identification", we draw 
three important conclusions: First, only the modulus of 
the form factors f{k) enters the noise correlations. Two 
types of order that differ only by a phase factor are thus 
difficult to distinguish by inspecting only the vicinity of 



the Fermi points. For instance, a bond-order- wave state, 
i.e., a bond-centered CDW with p-wave character, and 
a conventional site-centered 5-wave CDW lead to simi- 
lar features in the shot noise. Second, since the noise 
correlations are bounded by Eq. (4), they cannot di- 
verge. However, we expect any divergences present in 
the structure factors, properly rescaled, to also appear in 
the rescaled noise correlations LG, where L is the size of 
the system; see also Ref. 17. Third, it is interesting to 
note that only the nature of a fluctuation, i.e., whether 
it is particle- hole or particle-particle-like, determines the 
sign of the divergence. According to the analysis above, 
particle-hole correlations enter with a negative sign as 
dips along k' = k±q, while particle-particle fluctuations 
give a positive contribution and are observed along the 
anti-diagonal k' = —k of the shot noise. This is what 
is expected intuitively because a particle-hole scattering 
process is only effective if the /c-orbital is occupied in 
a different way than the k ± g'-orbital, i.e., if they are 
anticorrelated. On the other hand, for particle-particle 
scattering, the k and the — /c-orbitals must be both filled 
or both empty in order to be effective, i.e., they must be 
(positively) correlated. 



III. THE EXTENDEND HUBBARD MODEL 



We now consider the extended Hubbard model on 
a chain with L sites and periodic boundary conditions 
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(a) G^^(k,k'),L=50,N=50 (b) G^^(k,k') (c) G(k,k')/2 
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FIG. 1: (Color online). Noise correlations (a, b, c, e, f) and momentum distribution (d) in the extended Hubbard chain with 
IJ jt — V/t = 10 at filling x = 1/2, i.e., kp = 7r/2^ obtained using coordinate- space DMRG calculations and subsequent Fourier 
transformation. The system exhibits an ordered 2kF-CDW, manifesting itself in the pronounced dip along k' = k ±2kF in G^^ 
(a) and vanishing correlations in the ^[-channel (b). This can be best seen in the intensity plot of the total noise correlations 
(e). The rescaled shot noise (f) at k = kp, k' = —kp scales linearly with the size of the system L, indicating true long-range 
order. Note that because of the closed shell configurations used in these calculations, kp is not part of the discrete Brillouin 
zone. The finite-size scaling (f) thus shows paths along k' = —k + 2n/L. 



(PBC). The Hamiltonian of this system reads 

nini^i , (10) 

l,cr I 

where cj ^ (c^ ^) creates (destroys) an electron with spin 

a on site ni^a = a^i a occupation number 

operator and the parameters t', U and V characterize 
the next-nearest-neighbor hopping, the on-site and the 
nearest-neighbor interactions, respectively. We focus on 
the ground state sector with = = N/2, with N 
the total number of particles. The filling x is defined as 
X = N/(2L) and thus, as usual, half-filling {x = 1/2) 
corresponds to one electron per site on average. For con- 
venience, we set the hopping amplitude t = 1 and ex- 
press the other interactions in units of t. The extended 
Hubbard model is 5'/7(2)-invariant for all values of the 
couplings, and the ground state is a spin singlet for the 
parameter values considered in this work. Longitudinal 
and transverse spin correlations, which appear in G^^ 
and G||, respectively, are therefore identical by symme- 
try. Furthermore, the continuous U{1) gauge symmetry 



and the SU{2) spin invariance forbid the appearance of 
true superconducting or spin-ordered states in this ID 
model. However, the discrete lattice symmetries or the 
time-reversal symmetry can, in principle, be broken, giv- 
ing rise to commensurate charge-density or bond-order 
waves. 

For the generic case of repulsive interactions and in- 
commensurate fillings, the low-energy sector of the ex- 
tended Hubbard model is accurately described as a TL 
liquid [32] with four independent parameters that com- 
pletely specify the model: the velocities of the charge and 
spin excitations, Vp and v^r, respectively, and the correla- 
tion exponents Kp and K^j. Since the SU (2) symmetry of 
the extended Hubbard model implies K^j = 1, the asymp- 
totic behavior of the correlation functions depends only 
on Kp, which in terms depends on the interactions and 
the filling. In the absence of interactions, these veloci- 
ties reduce to the Fermi velocity vf^ and the correlation 
exponents are equal to 1. 

The Hamiltonian (10) has been studied extensively in 
the past (see, e.g., Ref. 22 for an overview). Based on this 
knowledge, we have chosen a few interesting parameter 
sets in the global phase diagram [23, 24, 25, 26, 27, 28, 
29] to investigate the noise correlations, see Tab. I. The 
presentation in Sec. IV is ordered according to increasing 
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FIG. 2: (Color online). Noise correlations (a, b, c, e, f) and momentum distribution (d) for the Hubbard chain with U/t — 10 
and V — at half -filling , i.e., kp — 7r/2. The system is in the Heisenberg regime where only spin-spin correlations are critical 
because the charge sector is gapped. Note that in (a), we have chosen a very small upper cutoff to reveal the structure around 
opposite Fermi points. The behavior along the diagonal can be deduced from (c). The system exhibits a 2kF-SDW, which 
manifests itself in both the tt (a) and the ^ I- channel (b) as dips along k' = k ±2kF- This is well visible in the intensity plot 
(e). In contrast to the case of Fig. 1, this system only exhibits quasi long-range order. The rescaled shot noise (f) thus grows 
much more slowly than in the case of true long-range order. 



complexity of the noise spectrum. 



IV. NUMERICAL RESULTS 

In this section, we present our numerical results for 
the noise correlations G (2) in the extended Hubbard 
chain (10) obtained by DMRG calculations in coordi- 
nate space and exact diagonalizations (ED) in momen- 
tum space [30]. Due to the vast number [0(L^)] of four- 
point correlation functions that have to be evaluated in a 
coordinate space approach, we only consider chains with 
up to 56 sites. These sizes are, however, sufficient for the 
purpose of the present work because the main features 
of the noise correlations are already well-established in 
these systems and because the finite-size scaling does not 
show any ambiguities. We have checked the implemen- 
tation and the accuracy of the coordinate-space DMRG 
calculations against ED results on smaller samples (up 
to 1/ = 20 sites) and have found very good agreement. In 
addition, our calculations all comply very well with the 
sum rule (5). 

The number of fermions N and the system size L are 
always chosen so that closed shelly configurations occur, 
i.e., so that all orbitals below kp are occupied, while or- 



bitals above kp are empty. This condition translates into 
a fermion number N = 2 + 4m, where m is a nonnegative 
integer and hence L = N/ (2x). On the one hand, this 
choice eliminates spurious effects arising from open-shell 
conditions, but, on the other hand, the Fermi vector kp 
cannot be part of the discrete Brillouin zone; only even 
multiples of it can be resolved. Nevertheless, the effects 
we will study will be well-visible in the vicinity of the 
Fermi momentum; its behavior for k ^ ^kp can easily 
be deduced. 



A. Ordered charge-density wave 

Let us start with a simple case that exhibits true long- 
range order by choosing a half- filled system {x = 1/2) 
with U/t = 10 and V/t = 10. In this case, the charge 
and the spin sectors are gapped. In the limit t = 0, 
V/2 > is a sufficient condition for the existence of 
a charge- density- wave state in which every second site is 
doubly occupied, see Tab. I. In this strong-coupling limit, 
the noise correlations can be calculated analytically using 
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FIG. 3: (Color online). Noise correlations (a, h, c, d, e) and momentum distribution (f) in the Hubbard chain with U/t = 10^ 
t' /t = 0.7; and V = at filling x = 1/2. The system is in a phase with long-range modulations of the kinetic energy. In 
G^^, (a), a rather sharp dip along k' = k ± tv can be seen. The finite-size scaling (d) illustrates the long-range character of 
the correlations. Note that because of the closed shell configurations used in these calculations, kp is not part of the discrete 
Brillouin zone. The finite-size scaling (d) thus shows paths along k' = —k + 2ti/L. In the ^[-channel (b), only a shallow 
valley can be identified along the diagonal k' — k^n. Finite-size scaling (not shown) indicates that these spin correlations are 
short-range. Due to the phase insensitivity, the noise correlations alone cannot distinguish between this bond-centered and a 
conventional site- centered CDW, when focussing only on a close neighborhood of the Fermi points. 



a site-factorized wave function 
Gtt (fc, fcV fc) = E E 

1 = 1 l' = l 

Gn{k,k') = 0. (11) 

Here ji is the period of the modulation and m an integer. 
We thus expect a signature along k' = k ± 27r//i. 

These arguments are in perfect quantitative agreement 
with the numerical results for the shot noise shown in 
Fig. 1, which clearly reveal a = tt CDW signature 
in the ||-channel and vanishing correlations in the ||- 
channel. Although the featureless momentum distribu- 
tion shown in Fig. 1(e) is an indicator of a gapped state, 
it does not reveal its nature. Since we have chosen the 
number of particles and the length of the chain to exclude 
the Fermi points, a slice along the anti-diagonal k' = —k 
does not contain the minimum of the dip. However, along 
a slightly shifted path along k' = — A: + 27r/L, as depicted 
in Fig. 1(f), it is easy to see that the minimum of the 
resettled shot noise increases linearly with the size of the 
system, in agreement with the strong-coupling derivation 



above. In the thermodynamic limit, the rescaled shot 
noise thus contains a diverging dip along k' = k ^ 2/ci?. 



B. Algebraically decaying spin-density wave 

Without the nearest-neighbor interaction V ^ one re- 
covers the standard Hubbard chain. At half-filling, the 
charge sector is gapped, but the spin excitations are gap- 
less. In this case, the correlation exponent Kp vanishes 
and density-density correlations decay exponentially. Be- 
cause of the SU (2) invariance, Kcr = 1 and the spin-spin 
correlations are critical, asymptotically decaying as 1/r. 
The low-energy sector of a half-filled Hubbard model is 
equivalent to a 5* = 1/2 Heisenberg chain. A snapshot 
of a spin configuration is depicted in Tab. I. The simple 
analytical strong-coupling analysis above is equally ap- 
plicable to spin modulations; the only difference is that 
spin correlations enter not only G^^, but also G||. The 
periodicity of the modulation is again two lattice spac- 
ings, leading to a signal along k' = k ±2kF- In contrast 
to the previous case, this system does not exhibit true 
long-range order because the SDW state would break a 
continuous symmetry. The ordering tendencies are ac- 
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cordingly less distinct than in the CDW example with 
true long-range order. This is illustrated in the finite-size 
scaling analysis of the rescaled shot noise LG in Fig. 2(f). 

In addition to the well-known staggered spin correla- 
tions, the Heisenberg chain also has staggered dimer- 
dimer correlations decaying with the same power law 
(disregarding logarithmic corrections). In fermionic lan- 
guage, these dimer correlations share the same properties 
as the bond-order- wave correlations, i.e., are modulations 
of the kinetic energy. In terms of higher angular momen- 
tum density waves, this is a p-wave CDW. We therefore 
see a superposition of both the spin and the kinetic en- 
ergy correlations in G^^ along k' = k ^ 2kF- 



C. Ordered bond-order wave (dimerized) state 

To extend the discussion on critical bond-order wave 
correlations in the half-filled Hubbard model, it is inter- 
esting to separate the two contributions by driving the 
system into a dimerized phase, where the kinetic energy 
is modulated and the spin-spin correlations decay expo- 
nentially. Such a phase has been reported in Ref. 29, 
in a study of the Hubbard model {V = 0) with an addi- 
tional next-nearest neighbor hopping t' ^ where a spin gap 
opens and the spin correlations become very short-range. 
In the large U limit, this model maps onto a frustrated 
Heisenberg chain, which is known to be dimerized at the 
Majumdar-Ghosh point J2/J1 = 1/2 [31]. In Ref. 29, it 
has been shown that this phase persists for finite values 
of U. We have therefore chosen to investigate the noise 
correlations for a half-filled chain with U/t = 10 and 
f /t = 0.7. In a state where the kinetic energy is modu- 
lated with period /i, we expect pronounced correlations 
of the form 



(4..c,+i,. + H.c.) (4 



.c.)) 



Tab. I depicts such a dimerized state (with modulation 
period /i = 2 lattice spacings) schematically. Writing 
the noise correlations (2) in coordinate space, one can 
convince oneself (after a bit of algebra) that a bond- 
order wave leaves a non-trivial fingerprint only in the 
||-channel along k' = /c±7r. Due to the particle-hole na- 
ture of the fiuctuation, the signature is negative. The ||- 
channel contains non-singular contributions of the form 
gizfcgizfc This crude approximation is confirmed by 
our numerical results presented in Fig. 3. The well- 
pronounced dip along /c' = /c ± tt in is clearly dis- 
tinguishable from the shallow surface in G||. The finite- 
size scaling [Fig. 3(d)] illustrates the long-range charac- 
ter of the correlations. In contrast, the finite-size scaling 
of G^i along the same path (not shown) quickly satu- 
rates, showing that the spin correlations are indeed short- 
range. Note that the noise correlations alone cannot dis- 
tinguish between this bond-centered and a conventional 
site-centered CDW in ID due to the phase insensitivity. 



D. Ordered charge-density wave coexisting with an 
algebraically decaying spin-density wave 

Another interesting example is a system with both 
SDW quasi-order and true long-range CDW ordering at 
different wave vectors, realized for a quarter-filled chain 
with U/t = V/t = 10. For these parameters, the model 
is an insulator with a charge gap, i.e., the correlation ex- 
ponent Kp = 0. As will be given explicitly in Eqs. (12), 
the spin-spin correlations decay as 1/r, as in the case 
of the Heisenberg chain. The strong-coupling ground 
state is shown in Tab. I. Based on this picture, we expect 
2kF spin modulations and Akp charge modulations to be 
present. This prediction is in excellent agreement with 
the numerical results of Fig. 4. A pronounced CDW re- 
sponse is visible along /c' = /c±7r in G^^ [Fig. 4(a)], while 
less pronounced SDW ordering can be identified in G^^ 
and G^i [Fig. 4(b)] along k' = k ± 7r/2. Note that the 
small dips along k' = /c ± 37r/2 in Fig. 4(b) are just due 
to the periodicity of the Brillouin zone. The finite-size 
scaling along the anti-diagonal clearly reveals that charge 
ordering is much stronger than the tendency to form a 
SDW, as can be seen in Fig. 4(f). 



E. Algebraically decaying singlet superconductivity 
and charge-density wave correlations 

For an attractive on-site interaction {U/t < 0), the 
fermions tend to pair in singlets, leading to a gap in the 
spin sector, i.e., it requires a finite energy to break up an 
on-site pair because, even in the thermodynamic limit, 
spin correlations and triplet superconducting correlations 
decay exponentially, while singlet superconducting and 
charge fiuctuations exhibit an asymptotic behavior of the 
form [23] 



Css (0 = (4^4^Co^^Co^^) ^ j^^^ , 

cos(2/cf0 



Ass 



CcDW (l) = (nino) ^ n ^Acdw- 



with amplitudes that depend on the values of the in- 
teractions and on the filling. The presence of a spin gap 
is taken into account by formally setting the correlation 
exponent Kcr 0. Since Kp > 1 for attractive inter- 
actions and away from half filling, the SS fluctuations 
are always dominant and the tendency to form a CDW 
is subdominant. Fig. 5 shows the noise correlations and 
the momentum distribution [Fig. 5(d)] calculated for a 
Hubbard model with attractive interaction U/t = —10 
(y = 0) at a filling of x = 3/8. Based on Ref. 22, we es- 
timate Kp ^1.3. Note that in Fig. 5(a), we have chosen 
a very small cutoff in order to display the much smaller 
CDW signature along k' = k^2kF in G||. As expected, 
SS correlations are visible along the anti-diagonal in G^i 
and clearly dominate over the CDW signature; this is evi- 
dent in the total shot noise. Fig. 5(c). In the region acces- 
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FIG. 4: (Color online). Noise correlations (a, b, c, e, f) and momentum distribution (d) in the extended Hubbard model 
with U/t — V/t — 10 at quarter- filling , i.e., kp — 7r/4. This system is characterized by long-range charge modulations at 
wave-vector Akp in G^^ (a) and algebraically decaying spin- spin fluctuations at 2kF in G^^ and G^i (b). The small dips along 
k' — /c ± 37r/2 are simply due to the periodicity of the Brillouin zone. A finite-size scaling analysis of the total shot noise along 
the anti- diagonal (f) nicely illustrates the difference between true and quasi-long-range order. 



sible to bosonization, i.e., around opposing Fermi points, 
our results for the noise correlations agree well with the 
behavior illustrated in the last column of Fig. 3 of Ref. 17. 
Note, however, that the sum rule (5) leads to nonsingu- 
lar structures in Such structures can be clearly seen 
in Fig. 5(f) as pronounced dips along k = -^kp- These 
dips and valleys are necessary to compensate the strong 
positive contributions due to the singlet superconducting 
correlations and are not predicted by bosonization. 



F. Metallic Tomonaga-Luttinger liquid 



We take the conventional repulsive Hubbard model at 
filling X = 3/8 with = 10 as an example of a spin- 
1/2 Luttinger liquid with both gapless charge and spin 
excitations. From Refs. 24, 34, we estimate Kp ^ 0.58, 
while = 1 due to the SU{2) symmetry. According 
to Refs. 25, 33, the density-density and spin-spin corre- 
lations decay with the same exponent 1 + Kp (neglecting 



logarithmic corrections) 



K, 



cos(2/ci?/) 



CcDW = (^^^o) n'^ — + ^CDW + l 



, cos(4/cfO 

^ CDW .,4^^ 



CsDW = {^i ^o) ^ ~ (^^ + ^SDW^^y^^ , (l^aj 

while the pairing correlations exhibit an unmodulated 
asymptotic behavior with an exponent 1 + 1/^p 



Css = (^4^c|^^co,|Co,t) ■ 



Bts 



(12b) 



The amplitudes A^.,B^ depend also on the interactions 
and on the filling. The shot noise obtained from this 
system is shown in Fig. 6. In contrast to the previous 
examples of insulating phases, the momentum distribu- 
tion [Fig. 6(d)] now has an infinite slope singularity at 
g = /ci?, as is characteristic of a TL liquid [35]. This pro- 
nounced feature can also be recognized in the noise corre- 
lations around the Fermi points, even for strong interac- 
tion strength. Although particle-hole fiuctuations (which 
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FIG. 5: (Color online). Noise correlations (a, b, c, e, f) and momentum distribution (d) in the attractive Hubbard chain with 
U/t — —10 (V — 0) at filling x = 3/8. In real space, electrons tend to form coherent on-site singlets. In momentum space, this 
pairing is reflected in dominant particle-particle correlations along k = —k' in G^i (b). In the total shot noise (f), it becomes 
apparent that the 2kF-CDW fluctuations present inG^^ (a) are decaying faster than the superconducting fluctuations . In contrast 
to the noise correlations, the momentum distribution (d) does not reveal any information on the ordering tendencies. Note in 
(f) the presence of a kink-like negative feature at -^kp, which compensates the positive response of the pairing correlations, by 
virtue of the sum rule (5). 



are always negative) are dominant for all cases of re- 
pulsive interactions, we expect pairing correlations (with 
positive contributions) to be visible as well, but with a 
weaker dependence on the system size. Our numerical 
results confirm these expectations: In Fig. 6(a) and 
G||, Fig. 6(b), the dips along k' = /c ± are well- 
pronounced, whereas the particle-particle peaks close to 
the Fermi points are barely discernible. It is interesting to 
note that the main features describing this metallic phase 
already emerge from a simple perturbative calculation; 
see appendix A and compare with Fig 7. The behav- 
ior around opposing Fermi points {k = kp^ k' = —kp) is 
in perfect qualitative agreement with recent bosonization 
results [Eq. (8) of Ref. 17]. The detection of Akp charge 
fluctuations, which are present in principle, is hindered 
by the fact that they decay as a power law with an expo- 
nent significantly larger than two and presumably have 
a small amplitude A'^j^^. 



G. Summary 

Let us now summarize the main signatures of different 
phases that can be observed in the noise correlations. 
As a general rule, particle-particle fluctuations give a 



TABLE II: Diverging behavior observed in the rescaled noise 
correlations LG for different phases. A plus (+) or minus (-) 
symbol indicates the sign of the singularity. 



Phase 




LG^l 


Charge-density wave (CDW) 







Spin-density wave (SDW) 






Singlet superconductivity (SS) 





+ 


Triplet Superconductivity (TS) 


+ 


+ 


Bond-order wave (BOW) / dimerization 








positive contribution and are observed along the anti- 
diagonal k' = —koi the shot noise, while particle- hole 
correlations enter with a negative sign and lead to dips 
along k' = k ± q. The strong signal along the diag- 
onal k' = k contains only information that is already 
encoded in the momentum distribution. Tab. II summa- 
rizes the five different phases encountered in the examples 
presented above, indicating the sign of the singularities 
in the rescaled noise correlations LG, together with the 
channels in which they can be observed. If there is only 
one ordering tendency in the system, it can be identi- 
fied unambiguously. However, as soon as there are two 
or more competing phases, one has to include additional 
information to distinguish between them. 
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FIG. 6: (Color online). Noise correlations (a, h, c, e, f) and momentum distribution (d) in the Hubbard model in a metallic 
phase with U /t =10 (V — 0) at filling x = 3/8. The metallic character is reflected by the infinite slope (not jump) singularity 
of the momentum distribution (d) at k = ±kF- This singularity also translates itself to the noise correlations. CDW and SDW 
signatures (with a negative sign, due to their particle-hole nature) are recognizable along k' = k±2kF, while the small positive 
peaks around the Fermi points are due to the pairing fluctuations. The chosen parameters correspond to a Luttinger liquid 
correlation exponent Kp ^ 0.58. 



V. EXPERIMENTAL CONSIDERATIONS 

Finally, we would like to put these numerical results in 
context with respect to experiments. Initially, the atoms 
would be confined in an optical lattice. By using PBC in 
our numerical calculations, we have neglected the pres- 
ence of a shallow confining potential as well as the open 
ends of the ID system. We have tested the sensitivity 
of the shot noise to the choice of boundary conditions in 
the case of a combined CDW/SDW; see Sec. IV D. Qual- 
itatively, we find exactly the same features for OPC as 
for PBC, namely, a pronounced dip along k' = k^Akp 
in G||, indicating the presence of a CDW, and a less 
distinct valley along k' = k ^ 2kF in G^^ and G||, re- 
fiecting the SDW ordering (not shown). After the trap 
is turned off, the atom cloud expands freely, allowing the 
momentum distribution n^r {k) as well as the noise cor- 
relations Gaa' {k^k') to be measured. The relationship 
between the freely expanding atom cloud and the initial 
lattice states of the trapped atoms is explained in detail 
in Ref. 11. We are thus confident that our results directly 
apply to experimental realizations and measurements. 
The summary of features present in the noise correla- 
tions (Tab. II) clearly illustrates the necessity for state- 
selective measurements to distinguish different phases, 
as the sign of the signal reveals only the character of the 



fiuctuation, i.e., particle- hole- or particle-particle-like. 



VI. CONCLUSION 

We have analyzed the noise correlations (density- 
density correlations in momentum space) in various 
phases of the extended Hubbard model in one dimen- 
sion. Our numerical density-matrix renormalization 
group study, carried out for several characteristic in- 
teractions and fillings, shows that different types of 
(quasi-) long-range order leave different fingerprints in the 
shot noise. This allows different phases to be identified 
and distinguished with a universal probe. The method 
is therefore of interest for experiments with ultracold 
atoms, in which the shot noise can be extracted from 
time-of-fiight images. We have also pointed out the im- 
portance of sum rules, i.e., the fact that the integral of 
over the Brillouin zone is equal to zero, which leads to 
the interesting effect that noise features due to correla- 
tions must be compensated by a (possibly nonsingular) 
complementary structure. In the future, it would be de- 
sirable to obtain similar ab-initio results for the noise cor- 
relations of various fermionic systems in two dimensions, 
notably the repulsive Hubbard model at and away from 
half filling. Furthermore the effect of finite temperature 
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(a) G^^(Kk')/lP 




(b) G^^(k,r)/U' 




(c) G(k,k')/(2LP) 




FIG. 7: (Color online). Noise correlations G/U'^ obtained in second- order perturbation theory for a metallic system at a filling 
of X = 3/8. When compared to the DMRG results shown in Fig. 6, the positive peaks in the vicinity of opposite Fermi points 
due to particle-particle correlations are more pronounced. The analytical expressions have been evaluated here for L = 56 and 
N 42. 



on the noise profiles needs to be investigated. 



Momentum distribution 
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At zero temperature, the momentum distribution of 
the noninteracting system reads 

n^^Hk) = f{k) = 9(\k\-kF) , 

where (k) is the unit step function. For convenience, we 
also define f {k) = 1 — f {k). The first corrections arise 
in second order perturbation theory and are therefore 
independent of the sign of U 



APPENDIX A: PERTURBATIVE REGIME 



n(2) (k) ^ n^^^ (k) + [h (k) - h {k)] 



(A4) 



For t' = V = ^ and small on-site interactions \U\ <C 1, 
one can calculate the noise correlations in perturbation 
theory. In the momentum representation, the Hamilto- 
nian (10) reads 



k,a 



k,k' ,q 



with the dispersion 

e/c = —2 cos /c + 2 coskp 



(A2) 



with 



h (fc) = 



h{k) 



fik)^ f{k')f{q)f{k + k' -q) 



E 



fe',q + ~ ^k+k'-q) 
q^k 

fik)^ f(k')fiq)fik'+q-k) 



E 



j^' q {^k' -^^q- ^k- ^k'+q-kY 



and Fermi momentum kp = xtt. To first order, the 
ground state wave function is given by 



I*) = 1^) + 7 E 



Ac, Ac ,g 

ki^k' 



^k^k' 



\F) , (A3) 



where \F) = Ylk<kF ^k a\^) Fermi sea. 



2. Noise correlations 

For non-interacting fermions on a chain with PBC, 
Gaa' {k, k') = because Uk^a commutes with the tight- 
binding Hamiltonian. Analogous to the momentum dis- 
tribution, the lowest order corrections to the noise corre- 
lations are quadratic in U and thus also independent of 
the sign 
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Sa,a'h,k'Il (k) + (1 - <5^,a') -^4 (fc, fc') 
-/3(fc,fc')+/l(fc)/2(fc') 
Sa,a'Sk,k'l2{k) + {l-6,,,>)l5{k,k') 



h (fc) /i (fc') {k>kFAk'> kp) 
{k > kp A k' < kp) 
h (k) h (fc') {k<kFAk' < kp) 



(A5) 



with the two-point correlations Ij given by The normahzed noise correlations GjU^ obtained in 

this perturbation calculation are shown in Figs. 7 for a 

I (j^ y\ _ f {q) f {k ^ q - k') metallic system at filling x = 3/8. Compared to the 

^ ^ ^ ^ 1^2 (^^^ ^ _ _ ^,^2 ' DMRG results shown in Fig. 6, the positive peaks in the 

^ vicinity of opposite Fermi points due to particle-particle 

I^ik k') — ^ ^ ^^'^ / {q) f {k ^ k' — q) correlations are more pronounced, but, otherwise, the 

' ^ (e/e + ek' — — ^k+k'-qf ' pcrturbativc approach leads to a surprisingly accurate 

, , , _ ^ _ ^ , ^ description of the shot noise. 
T f{q)f{k^k'-q) 

^5 Kf^.i^ ) - Y2 2^ - ^ • 

Q \^k + ^k+k'-q — ^k — ^k' ) 
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